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An arbitrarily reliable quantum computer can be efficiently constructed from noisy com- 
ponents using a recursive simulation procedure, provided that those components fail with 
probability less than the fault-tolerance threshold. Recent estimates of the threshold are 
near some experimentally achieved gate fidelities. However, the landscape of threshold 
estimates includes pseudothresholds, threshold estimates based on a subset of compo- 
nents and a low level of the recursion. In this paper, we observe that pseudothresholds are 
a generic phenomenon in fault-tolerant computation. We define pseudothresholds and 
present classical and quantum fault-tolerant circuits exhibiting pseudothresholds that 
differ by a factor of 4 from fault-tolerance thresholds for typical relationships between 
component failure rates. We develop tools for visualizing how reliability is influenced by 
recursive simulation in order to determine the asymptotic threshold. Finally, we conjec- 
ture that refinements of these methods may establish upper bounds on the fault-tolerance 
threshold for particular codes and noise models. 
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1 Introduction 

A quantum computer can potentially solve certain problems more efficiently than a clas- 
sical computer m^TEl- However, quantum computers are likely to be engineered from inher- 
ently noisy components, so any scalable quantum computer system will require quantum error 
correction and fault-tolerant methods of computation. As candidate quantum device tech- 
nologies mature, we need to determine component failure probabilities necessary to achieve 
scalability. The fault-tolerance threshold for gate and memory components is particularly in- 
teresting because arbitrarily reliable computations are possible if the circuit components have 
failure rates below the threshold. Given detailed knowledge of the fault-tolerance thresh- 
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old and its associated trade-offs, proposals for fault-tolerant quantum computation can be 
critically evaluated. 

The concept of a fault-tolerance threshold has its origins in the classical theory of compu- 
tation. In the 1950's, von Neumann showed that it is possible to achieve a reliable classical 
computation with faulty components provided that the failure probability of each component 
is below some constant threshold probability that is independent of the circuit size and the 
desired noise rate 0j. Similarly, concatenated coding and recursive error correction can be 
used to achieve reliable quantum computation. Concatenation is the process of encoding 
physical bits of one code as logical bits of another code. It is now well-known that using a 
single-error-correcting concatenated coding scheme with L levels of recursion, the maximum 
failure probability "/circuit of a fault-tolerant circuit can be estimated as a function of the max- 
imum failure probability 7 of a basic component using the fault-tolerance threshold inequality 



Jcircuitil) < f 7 V 

7th ~ \lth) 



(1) 



where jth is the asymptotic threshold [S]. When Eq JIJ holds with equality, we call it the 
fault-tolerance threshold equation. The final circuit failure probability "/circuit decreases as a 
doubly exponential function of L if 7 < j t h- 
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Fig. 1. Plots of quantum threshold estimates for stochastic noise models between the years 1996 
and 2004 E]|S|[7||HliniEniE]ElElC3El- Stars denote numerical estimates and circles denote 
all others. Most of these estimates apply to the [[7, 1, 3]] code, though the pair of thresholds at 
and above 10~ 2 apply to surface code memories and post-selected quantum computation. The 
networks vary over unitary networks, nearest-neighbor networks, and various optimized networks. 
Some estimates apply only when there are no memory errors, and others apply only for the Clifford- 
group gates. The dark swath designates the large interval that contains all of these estimates. 

One branch of fault-tolerant quantum computing research has focused on estimating the 
fault-tolerance threshold in the fault-tolerance threshold in equality . Figure shows the range 
of quantum threshold estimates reported in [SI El E| El El HD1 ID1 C2l 1131 DHI CSI - The estimates, 
which vary between 10 -6 and 10~ 2 , include numerical and analytical results for varying 
networks for the [[7, 1, 3]] code, such as optimized networks, unitary networks, and nearest- 
neighbor networks. They also include results for surface code memories and post-selected 
quantum computing models that yield thresholds above 10 -2 . Numerical estimates tend to 
be more optimistic than their analytical counterparts. 
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Fig. 2. (a) An ideal threshold reliability information plot (TRIP) follows Eq The crossing 
point between the L = line and the L = 1 curve, marked by the asterisk on the thick vertical 
line, is the fault-tolerance threshold. All of the curves cross at the same point, (b) A real TRIP 
does not follow Eq 0. The crossing points between the L = line and the other curves are 
all different. These points, marked by circles, correspond to a sequence of pseudothresholds that 
converges to the real fault-tolerance threshold marked by the asterisk on the thick vertical line. 



If there is a single type of component that is replaced using the same rule each time, 
the fault-tolerance threshold inequ ality, Eq QJ, becomes an equality. This leads to a very 



simple prediction shown in Figure 2(a) if 7 = 7^, then "/circuit — 1th- This fact is taken 
as the basis of several numerical analyses of the fault-tolerance threshold today. Specifically, 
this simplification is attractive for computationally expensive numerical simulations because 
it implies that the threshold can be determined by finding the smallest nonzero value of 7 
that solves ^drcuitil) = 7- 

Realistically, however, there are multiple types of components that are each replaced using 
different rules, so the f irst crossing point does not accurately indicate the fault-tolerance 
threshold. Figure 2(b) more accurately portrays the effect of recursive simulation. As the 
recursion level increases, for example, an exponentially growing number of wires must be 
introduced between gates. When these wires are unreliable, as they likely will be in quantum 
circuits, successive recursion levels can cause errors to increase even though 7 is beneath the 
apparent threshold. Thus, recursive simulation changes the relative proportions of each type 
of component and what appears to be the threshold at one level of recursion may be far 
from the asymptotic threshold. This sequence of crossing points cannot be used to describe 
the proper conditions under which a system is scalable. Rather, these crossing points are 
pseudothresholds (Hj- 

When multiple types of components are replaced using different rules, each component 
type must be parameterized by a separate failure probability. Hence, opportunity exists for 
engineering trade-offs that still preserve scalability. These trade-offs can be quantified given 
the asymptotic threshold, families of pseudothresholds, and their relation to the shape of the 
set of subthreshold component parameters. 

In this paper, we present practical methods for distinguishing pseudothresholds from 
asymptotic thresholds. In particular, we explore the conditions under which pseudothresholds 
exist and clarify their meaning. We em bark on this exploration carrying two tools. The first 
tool is a generalization of Figure [2(b)| which we call the threshold reliability information plot 
(TRIP). In a TRIP, each curve represents the failure probability of a particular component 
at concatenation level L and crosses the L = line once. The crossing of a level-!/ curve 
and the lcvcl-(i = 0) line yields the rightmost edge of an interval on the 7-axis below which 
reliability is improved by concatenation. The crossing point is a level-L pseudothreshold. 

The second tool is a threshold information flow diagram (TIFD) that shows how recursive 
simulation can change the reliability of a particular set of noisy components. A flow is 
a normalized vector field that can be visualized as a collection of arrows. Each arrow's 
base is anchored to a point that represents the current failure probability of the (composite) 
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Fig. 3. Threshold information flow diagram (TIFD) corresponding to the following recursive 
simulation procedure: gates u and v at level-(L — 1) are both replaced by fault-tolerant gates at 
level L that can withstand a single level- (L — 1) gate failure. However, u and v compute different 
functions, so their fault-tolerant implementations are different. In this example, the fault-tolerant 
u contains two u gates and two v gates, and the fault-tolerant v contains three u gates and three v 
gates. Arrows on this TIFD indicate how the recursive simulation procedure changes the effective 
failure probabilities of u and v. For some failure probabilities, recursive simulation has no effect. 
These fixed points are marked by circles, and one of them determines the fault-tolerance threshold, 
indicated by a thick black line. Sample trajectories begin at the squares and flow along the thick 
dashed lines, where the diamonds mark the sequence of points for each level L of concatenation. 
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components. The direction of each arrow indicates the direction the anchor point moves 
at the next level of recursion. In contrast to the TRIP, the TIFD exposes how all of the 
component failure probabilities change in a recursive simulation. If the recursive simulation 
is self-similar, i.e., if the failure probability of a level- L component can be expressed in terms 
of the failure probabilities of level- (L — 1) components, particularly with respect to the noise 
model, then the TIFD indicates whether or not recursive simulation increases or decreases 
each component's failure probability, allowing us to compute and visualize the fault-tolerance 
threshold. 

For example, Figure shows a TIFD for a hypothetical pair of faulty basic gates u and v. 
Because u and v compute different functions, their fault-tolerant implementations use different 
numbers of basic u and v gates. In this example, the fault-tolerant u contains two u gates and 
two v gates, connected in some fashion so that the fault-tolerant u can withstand one internal 
gate failure and still produce a "good" output. If basic gates u and v fail independently with 
probabilities 7„ and 7„, respectively, then the failure probability of the fault-tolerant u gate 

l-(l-7«) 2 (l-7«) 2 -27n(l-7«)(l-7.) 2 -27,(l-7 t -)(l-7u) 2 - (2) 

Similarly, the fault-tolerant v gate contains three u gates and three v gates and can withstand 
any single failure, giving 

1 - (1 - 7u ) 3 (l - 7,) 3 " 3 7 «(1 - 7n) 2 (l - 7,) 3 " 37,(1 - 7«) 2 (1 - Inf. (3) 

There are many such hypothetical examples. In Section we give a more realistic example 
that demonstrates how a TIFD is calculated based on an actual circuit. 

Continuing with this example, the effective failure probabilities of u and v after level- 
1 recursive simulation both depend on the initial failure probabilities 7 ° and 7° of u and 
v, shown on the horizontal and vertical axes of the TIFD. Consider the following scenario. 
The v gate initially fails with probability 0.2 and the u gate does not fail at all. A square 
marks this point on the TIFD. The arrow at this point on the TIFD points down and to the 
right, indicating that a level- 1 recursive simulation will improve v but make u worse. The 
dashed line connects the initial failure probabilities to the failure probabilities of the level- 

1 simulated gates (at about (0.04,0.1)). The dashed path shows that subsequent recursive 
simulation makes u and v arbitrarily reliable. 

The TIFD in Figure El also indicates the set of initial failure probabilities that is below 
threshold. The boundary of this set is determined by a fixed point of the recursive simulation 
procedure. This fixed point is marked with a circle, and the thick dark line passing through 
this circle is the invariant set that indicates the fault-tolerance threshold. For example, an 
ideal v gate and a u gate that initially fails with probability 0.28 (marked by a square) is 
above threshold. This point flows nearly parallel to the invariant line at first, but ultimately 
escapes from the boundary of the TIFD after about 7 levels of recursive simulation. 

We organize the paper as follows. In Section[21 we first define the fault-tolerance threshold 
for concatenated flow maps, then we define pseudothresholds and describe their importance. 
We calculate, in Section a family of pseudothresholds and distinguish them from the fault- 
tolerance threshold estimate for the classical repetition code. In Section we study the 
[[7,1,3]] CSS code, again comparing the fault-tolerance threshold against families of pseu- 
dothresholds. We suggest techniques for finding the threshold in Section [5] that expand upon 
our use of the TIFD. We conclude in Section HJ with open questions. 

2 Pseudothresholds 

In Section^ we discussed the phenomenon of pseudothresholds without introducing many 
mathematical concepts. In this section, we clarify what we mean by defining both the fault- 
tolerance threshold and the pseudothreshold sequences for a set of flow maps. 

2.1 The fault-tolerance threshold 

The threshold is related to the number of "bad" fault paths through a circuit. Assume 
faults occur at discrete locations and times with probability dependent on the location type. 
The failure probability of the circuit can then be determined as a function of failure proba- 
bilities of types of locations. For a code correcting t errors, the probability that one of these 
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bad fault paths occurs is no more than Cj t+1 , where 7 is the largest failure probability of 
any type of location in the circuit and C is the number of ways to choose t+1 failed locations 
out of N total locations. The fault-tolerance threshold satisfies 

1th > f , (4) 

since ^th is the fixed point of the map "fdrcuitij) = C'j t+1 . 

However, we must recognize that the fault-tolerant implementations of each location type 
differ, so we must express the failure probability of each fault-tolerant location type as a 
function of the location types it contains. In other words, we construct (approximations 
to) the flow maps for the given fault-tolerant implementations and noise model jHJ 12 ■ In 
particular, if each type of location I is assigned an initial failure probability 7° and if there 
are n different types of locations, the approximate failure probability of location type £ after 
one level of recursive simulation is a function T\ of all n of the initial failure probabilities. 
Therefore, T\ is called the flow map for location type I. 

Considering all types of locations I, the functions T\ are the coordinates of a flow map T , 
The flow map takes the failure probabilities of the n location types to their new values after 
one level of recursive simulation. The failure probabilities of the L-simulated location types 
are (approximately) related to the initial failure probabilities 7° by the composed flow map, 

r L » r 1 o. . -or' . (5) 

L times 

Ideally, the replacements can be constructed so that Eq (J5J is an equality This is the case for 
the example in Section |21 but not in Section 0] Let Tf denote the coordinate function of T L 
associated with location £, and let T° be the initial function that selects the I coordinate. In 
other words, Tf is the failure probability of location I after L levels of recursive simulation. 
The function Tf is a concatenated flow map for location type I because we can use Eq (JjJJ to 
derive (an approximation to) Tf. The map T L is the concatenated flow map. 

A vector of failure probabilities 7 for the n location types is below threshold if all n of the 
failure probabilities approach zero as the concatenation level approaches infinity, 

hmr L (7)=0. (6) 

L — >oo 

Let T be the set of these vectors that are below threshold, and let C c be the n-dimensional 
cube of edge length e with one corner at the origin, 

C e ee {7 e [0, 1]" \^<e W,e 6(0,1]} (7) 

The cube contains all of the vectors whose worst failure probability is less than e. The fault- 
tolerance threshold or asymptotic threshold is the size of the largest cube contained in T, i.e., 

lth ee sup{e > I C t C T}. (8) 

If all component failure probabilities are beneath this probability, then composing the flow 
maps reduces the failure probability arbitrarily close to zero. 

2. 2 Definition of pseudothresholds 

Before we define pseudothresholds, we introduce the concept of a setting. Settings param- 
eterize a set of location failure probabilities by a single parameter so that we can think of 
Tf as a function of this parameter. A setting is a function from a single failure probability 
parameter to a vector of n failure probabilities, one for each location. For example, the di- 
agonal setting 5(7) = (7, . . . , 7) is used in analyses that assign each location the same initial 
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failure probability. The Steane setting 5(7) is another setting that se ts all location failure 
probabilities to 7 except for a waiting bit which is assigned 7/IO ( m this reference, a 
waiting bit is assigned a failure probability of 7/ 100). 

Suppose there are n types of locations. We define a pseudothreshold 7^ for a fixed level 
of recursion L > 0, a location £, and a setting g as the least nonzero solution to 



17(5(7)) =7, (9) 

This definition presents the (L, i, g)-pseudothreshold as a fixed-point calculation for a function 
derived from the flow m ap T. The left-hand size of Eq @ can be viewed as one of the curves 



plotted in Figure 2(a) and the right hand size can be viewed as the L = line. The point 
where these curves intersect is a pseudothreshold. 

For fixed location £ and setting g, the sequence is not necessarily constant as a function 
of L. In fact, the sequence is typically not constant, meaning that pseudothresholds are a 
generic phenomenon. More specifically, let 70 be any pseudothreshold of the flow map T 
for a setting g, and let 70 be the constant vector (70, . . . , 7o)- The pseudothreshold 70 is 
independent of location type and recursion level only if the setting satisfies r (17(70)) = 70 and 
70 is a fixed point of T. 

Despite the fact that pseudothresholds are not thresholds, pseudothresholds are interesting 
because only a fixed level of recursive simulation will be used in practice. If the goal is to 
construct a reliable fault-tolerant location type £, and all of the location types have the 
same initial reliability (i.e., 5(7) = (7, . . . ,7)), then choosing 7 to be less than the (L,£,g)- 
pseudothreshold makes the L-simulated gate location type £ more reliable than the initial 
gate. However, some caution must be applied to pseudothresholds as well because the (1, £, g)- 
pseudothreshold and the (2, £, (^-pseudothreshold can be substantially different. 

In the following sections, we present an illustrative example of classical pseudothresholds 
followed by a more detailed example of quantum pseudothresholds. We show by means of 
these examples that pseudothresholds are generic to all multiparameter maps. In addition, 
we highlight that threshold estimates should account for multiple location types and higher 
levels of code concatenation to achieve more realistic threshold results. 



3 Classical Pseudothresholds for the [3,1,3] Code 

In this section, we analyze a classical example to build intuition about the differences be- 
tween pseudothresholds and thresholds. We study pseudothresholds for classical fault-tolerant 
components based on the [3, 1, 3] repetition code. We use the threshold reliability information 
plot (TRIP) of the [3, 1,3] code to identify pseudothresholds. We then characterize the flow 
map for this example using a threshold information flow diagram (TIFD). 

3.1 The [3, 1,3] code and its failure probability map 

In this example, the classical single-error-correcting [3, 1, 3] repetition code, also called 
triple modular rendundancy (TMR), is used to encode a single bit in three bits by copying- 
it three times. To make a fault-tolerant classical wire using this code, three location types 
Q = {u>, v, /} are required, where 

• w : {0, 1} — > {0, 1} defined by w(a) — a is a wire. 

• v : {0, l} 3 — > {0, 1} defined by v(a, b, c) = ab © be © ca is a voter. 

• / : {0, 1} — > {0, l} 3 defined by f(a) = (a, a, a) is a fanout. 

The superscripts above indicate the Cartesian product. The wire w is analogous to a wait- 
ing bit in a quantum fault-tolerance analysis, and the voter v and fanout / perform error 
correction. 

A noisy version of each location type is defined as follows. A noisy wire flips the output 
bit with probability j w . A noisy voter incorrectly indicates the output bit with probability 
7„. For simplicity, we choose to model the fanout gate to be noiseless. 
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To recursively construct a fault-tolerant wire, replacement rules are used. A replacement 
rule is a pair (b, R(b)) where b £ ft and R(b) is a circuit over that specifies how to replace 
a level-(L — 1) location at level L. R(b) is called the replacement of b and must preserve the 
functionality of the original location b. We construct the replacement rules to mirror replace- 
ment rules for quantum circuits, in which a circuit location is replaced by error correction 
followed by a fault-tolerant implementation of the location. 

The following steps suggest how to ensure proper component connectivity for a code 
encoding a single bit. First, replace b directly by D® n ° o R(b) o E® ni where D and E are an 
ideal decoder and encoder. We must have D o E equal to the identity gate on a single bit, 
where the open circle o denotes function composition. The numbers and n Q are the number 
of input and output bits of 6, respectively. After replacing each component in this manner, 
make a second pass over the circuit and replace all pairs D o E by bundles of wires. Finally, 
replace the remaining encoders and decoders by respective fault-tolerant implementations of 
input preparation and output readout. The resulting circuit components will be properly 
connected, and the circuit will not contain decoding and re-encoding components since these 
components are not typically fault-tolerant. 

For this example, a wire w is replaced by error correction followed by a transversal im- 
plementation of the wire, i.e., a wire is applied to each bit of the encoded input, shown in 
Figure^] Note the first dashed box indicates the classical error correction, which involves w, v, 
and / location types, and the second dash ed bo x in dicate s the fault-tolerant implementation 



of the original location. Similarly, Figures 5(a) and 5(b) show the fault-tolerant replacement 
of v and /, respectively. 




Fig. 4. Replacement rule for a wire w. The fanouts / and voters v perform error correction. 
The first dashed box indicates classical error correction using fanouts / and voters v. The second 
dashed box indicates the fault-tolerant implementation of the wire w. 



Using the replacements R{w) and R(v), the failure probabilities T 1 (7) and r*(7*) can be 
found, where the initial vector of failure probabilities is 7 = (7^,7°). Failure is defined to 
occur when the component's output does not decode to the correct value. The wire failure 
probability is easily calculated by counting the number of ways each configuration of errors 
occurs. For example, three voters fail in one way, three voters and one wire fail in three ways, 
two voters fail in three ways, etc. Rewriting the resulting polynomial in distributed form gives 

rl(7) = 67,7™ + 37" + 3 7 ^ - H - lfrtfy* + 127, 3 7™ ~ 187,7' 

+ 367^2 - 24 7 3 7 2 - 2 7 ,3 + I2 lvl l - 24 7 2 7 3 + IQ^l , (10) 

where the superscript has been dropped for notational convenience. 

Replacing the wire failure probability 7^ by 7" and the voter failure probability 7" by the 
probability of error correction failure 3(7°) 2 (1 — 7°) + (7°) 3 gives 

Ti(7) = H + 16 7 „ 3 - 39 7 „ 4 - 126 7 „ 5 + 474 7 „ 6 - 288 7 „ 7 - 936 7 „ 8 

+ 2080 7 ^ - 1824 7 i° + 768 7w n - 128 7w 12 , (11) 

where the superscript has been dropped again in the last expression. 
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Fig. 5. (a) Replacement rule for a voter v. The first dashed box indicates classical error correction 
using fanouts / and voters v. The second dashed box indicates the fault-tolerant implementation 
of the voter v. (b) Replacement rule for a fanout /. Because we assume fanouts are noiseless, the 
replacement is just three fanout gates /. 



The flow maps Eq IjlOl) and Eq (|llfl are exact and satisfy Eq JSJ) with equality, so they 
contain enough information to determine the asymptotic threshold. They can also be used 
to determine a bound on the number of "bad" fault paths, as discussed in Section [2 by 
considering only the low-order terms in the flow maps: 

ri,(7) < 3 7 ^ + 3 7 „ 2 + 6 7 »7. (12) 
Ti(7) < H + 167, 3 = (3 + 16 7 ,)7, 2 - (13) 

These bounds clarify the relative contribution each component makes to the failure probability 
of a fault-tolerant component. They also suggest a conservative bound of 1/12 on 7^. This 
can be calculated by assuming that 7 „ = "f w , solving for the fixed point of the right hand side 
of each inequality, and taking the least such fixed point . 

3.2 TRIPs for the [3, 1,3] code 

What is the behavior of the wire and voter failure probabilities as the concatenation level 
L increases? TRIPs based on Eqs (|1()|) - (|11|) provide a visualization of each level crossing point 
for the two types of locations. 

Figures 

ffl and □ are TRIPs for the voter and wire locations, respectively. From Figure El 
it is clear the voter probability r*(7) is a one-parameter map, and thus should resemble the 



TRIP for the fault-tolerance threshold equation (Figure 2(a) I. Since rj,( 7 ) is a function of 

only 7° for all L, each r^( 7 ) intersects at the same fixed-point of the map. This fixed point 
7th, indicated by a thick vertical line, is the L = 00 pseudothrcshold and the asymptotic 
threshold. It occurs at approximately 0.246. Note that the L — 00 pseudothrcshold for a 
particular location and the asymptotic threshold are not always equal, but they happen to 
be in this example. 

Even for a classical setting, there is a difference between pseudothresholds and the asymp- 
totic threshold jth- Figure Efshows the TRIP for the wire location at levels L = 0,l,2,3,oo. 
Here, unlike in the TRIP for the voter location, pseudothresholds appear in addition to an 
asymptotic threshold. This is because the replacement R(w) for a wire includes locations of 
type w, v, and /, creating a two-parameter map (fanouts are noiseless) that exhibits changing 
behavior with each concatenation level L. The curves now cross the L = line at differ- 
ent points. Each of these crossing points in the TRIP is a level-L pseudothreshold. As we 



A flow-map model for analyzing pseudothresholds 




0.1 0.2 0.3 0.4 0.5 



r° 

V 

Fig. 6. TRIP for a voter location for the [3, 1, 3] code for L = 0, 1, 2, 3, oo. Because is a function 
of only 7° constructed by recursive application of R(v), all of the curves intersect at the same 
point. This point is the fixed-point 7^ as 0.246 of the map and is indicated by an asterisk. 
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Fig. 7. TRIP for a wire location for the [3, 1, 3] code for L = 0, 1, 2, 3, 00 using the diagonal setting 
lw = 7?- Because is a function of both 7^ and 7JJ constructed by recursive composition, the 
curves cross the L = line at different points. Each of these points is a level- L pseudothreshold 
indicated by a circle, and the sequence of pseudothresholds 7^, 7^,, . ■ ■ , converges to the fault- 
tolerance threshold 7^ 0.246 indicated by an asterisk. 
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Fig. 8. TIFD for -f m and f v for the [3,1,3] code. The fixed points are indicated by circles. 
The region [0, 1/2) X [0,fth) enclosed by a thick line is the set of points below threshold, where 
7th ~ 0.246. Two sample trajectories begin at the squares and flow along the thick dashed lines, 
where the diamonds indicate the sequence of points as L increases, toward (0,0) and (1/2, 1/2). 



repeatedly replace the wire using R(w) 7 the number of voter locations begins to dominate, 
so the crossing point approaches the voter threshold, which corresponds to the asymptotic 
threshold. The difference between the asymptotic threshold jth ~ 0.246 and the level- 1 wire 
location pseudothreshold 7^ w 0.129 is 0.117, or 7^ w 1.9 x 7^. 

3.3 TIFDs for the [3, 1,3] code 

Given that pseudothresholds can be so different from 7^, can jth be determined from just 
one application of the flow map? In Section it was suggested that a TIFD provides an 
informative view of the effect of recursive simulation. Although a TRIP provides a visualiza- 
tion of the asymptotic behavior, it hides the fact that T 1 acts on a multidimensional space. 
When r 1 is chosen such that Eq (JSJ) is an equality, T 1 contains all of the information about 
the flow since T L can be expressed in terms of T 1 . A TIFD shows where each point flows 
under repeated application of T 1 without iterating the map explicitly. 

In Figure^ a TIFD for wire and voter failure probabilities on the unit half square is shown. 
The arrows represent the vectors T 1 (7) — 7, which give the probability flow under recursive 
simulation. There are five fixed points of the map: (0,0), (1/2,0), (1/2,7^), (1/2, 1/2), and 
(1,0), where 7^ w 0.246 cannot be expressed in radicals. Circles mark four of these fixed 
points. The subthreshold region T is [0,1/2) x [0,7th), indicated by the thick black box. 
Three corners of the subthreshold region are fixed points of the map. The fault-tolerance 
threshold j t h is the size of the largest "cube" , a square in this case, that is contained in the 
subthreshold region. 

First, if 7° = and 7° > 0, then the flow draws these points off of the 7°-axis. This 
occurs because the wire replacement rule R(w) contains both voter and wire locations, so the 
failure probabilities "mix". Second, if 7° = 0, then any point 7^, < 1/2 flows to the origin 
because the voters amplify any bias toward success. Third, the voter failure probability 
1^(7) is independent of 7^, so the voter probability is a simple one-parameter map under 
replacement. If 7° < jth, the map's fixed point, then the voter probability flows toward the 
7°-axis. Finally, if 7° < jth and 7^ < 1/2, then initially the wire failure probability may 
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increase because the voters are not reliably correcting errors. However, the voters improve 
with each iteration, so eventually this trend reverses. The voters begin correcting more errors 
than they introduce, and all of these points flow toward the origin. 

For the classical [3,1,3] code, the TIFD fully characterizes the threshold set T for three 
reasons. First, deviated inputs, i.e., inputs that are not codewords, are corrected before 
faults are introduced by the simulated gate locations. Second, the flow maps are the precise 
component failure probabilities. Third, there is no phase noise in classical fault-tolerance, so 
the parameters of the noise channel only change in the trivial way. Furthermore, the entire 
flow is easily visualized since the TIFD is two-dimensional. Under these conditions, the TIFD 
is an ideal tool for understanding and visualizing the process by which recursive simulation 
improves reliability and exhibits a threshold. 

4 Quantum Pseudothresholds 

We have seen that pseudothresholds exist even in a simple classical fault-tolerance scheme. 
In a quantum fault-tolerance scheme, the tools we have developed can now be applied to 
determine sequences of quantum pseudothresholds. In this section, we study thresholds for 
quantum fault-tolerance using the [[7,1,3]] CSS code. We follow the circuit construction 
given in ^2]. As in the classical example, TRIPs are again used to identify pseudothresholds 
for the given location types using particular settings, allowing us to determine the reliability 
achieved with each level of concatenation. In addition, we characterize the failure probability 
map using TIFDs. 

4-1 The [[7,1,3]] code and its flow map 

The [[7,1,3]] quantum code encodes a single qubit in 7 qubits with distance 3, meaning 
it corrects a general quantum error on a single qubit. The set of location types £ involved in 
the error correction routine is f2 = {1, 2, w, lm,p}: 

• £ = 1: one-qubit gate 

• £ = 2: two-qubit gate 

• £ = w: wait (memory) location 

• £ = lm: one-qubit gate followed by a measurement 8J, since the replacement rule for a 
measurement contains no error correction and since measurements in the networks for 
the 7-qubit code only appear after one-qubit gates 

• £ = p: ancilla preparation location, which we model as a one-qubit gate for simplicity 


We consider the depolarizing error model, where a location £ fails independently with 
probability 7^. In our probabilistic error model, we assume for a location failure on a single 
qubit, a X,Y, or Z error occurs with probability 7^/3. The distinction between X,Y, and Z 
errors is used in the approximation analysis. 

As in the classical case, we define a replacement rule for each type of location. We use 
the same replacement rules as given in , where we replace each location by error correction 
followed by a fault-tolerant implementation of the location. For example, a one-qubit gate is 
replaced by error correction and a transversal one-qubit gate, as shown in Figure El 

We use the same approach to derive the failure probability map as in the classical example, 
with the caveat that we do approximate counting of the ways in which two or more errors 
occur on the data qubits during the execution of the concatenated circuit. The composition 
of this map approximates the behavior of the concatenated circuit. This means that threshold 
results derived from this map are also approximate. The map does not account for incoming 
errors since when failure probabilities are below threshold, the probability of incoming errors 
should typically be small, and thus should not greatly affect the probability of two or more 
errors on the data qubits. The details of the failure probability map are given in Section IV 
Bof0. 
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Fig. 9. Replacement rule for a one-qubit gate U. The replacement includes error correction (EC) 
followed by a fault- tolerant implementation of the one-qubit gate U, indicated by a dashed box. 



4.2 TRIPs for the [[7,1,3]] code 

To determine pseudothresholds, we plot the reliability of each component at each level of 
concatenation using the flow maps. Figures HUH12I show TRIPs for wait, one-qubit, and two- 
qubit locations for the Steane setting (7(7) = (7,7, .. . ,7,7/10), where the last component is 
the wait location failure probability. 

These results in this section are obtained using a flow map, but they are verified using 
a Monte-Carlo simulation for one level of code concatenation. The Monte-Carlo simulation 
method randomly generates faults within the fault-tolerant circuit representing the original 
location according to the failure probabilities 7^. These faults create errors that propagate to 
the output of the circuit. A particular set of faults may or may not create too many errors at 
the circuit's output. By running on the order of one million simulations, we can estimate the 
failure probability of the original location. The particular Monte-Carlo simulation includes 
the ancilla preparation networks but does not model input errors so as to agree with the 
assumptio ns u nder which the flow map is derived. 

Figure ED snows the TRIP for a wait location at levels L = 0, 1, 2, 3, 00. Since the flow map 
is a multi-parameter map, the crossing points no longer cross the line L = at the same point 
and thus pseudothresholds appear at each level. As L increases, the level- L pseudothreshold 
approaches an asymptotic threshold that depends on the location and setting. Note that the 
region between the L — 1 curve and the L — line is quite small for these initial conditions, 
smaller than the region considered to be below the asymptotic threshold. This is similar 
to the behavior of the classical wait location. The behavior is largely due to the Steane 
setting - the level- 1 simulation of the wait location includes other location types that have 
been set to fail with probability 10 x 7^. However, as the level of concatenation of the wait 
location increases, the trade-offs between the failure probabilities of the location types begins 
to stabilize causing the level-4 pseudothreshold, for example, to be much larger and closer 
to the asymptotic threshold. The level-1 pseudothreshold as calculated from the flow map 
is 2.2 x 10 -5 and the corresponding value from Monte-Carlo simulation is approximately 
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Fig. 10. TRIP for a wait location for L = 0, 1, 2, 3, oo for the initial setting 7JJ, = 1/10(7?), where 
i indicates all location types except for a wait location. Circles indicate pseudothresholds and an 
asterisk marks the threshold for this gate and setting. The pseudothresholds occur at probabilities 
2.2 X 10~ 5 , 3.0 X 10~ 5 , and 3.4 X 10~ 5 . The wait threshold for this setting is 3.6 X 10~ 5 . This 
corresponds to 7 = 3.6 X 10~ 4 in the Steane setting. 



2.4 x 1CT 5 . 

Figure ITTI shows the TRIP for the one-qubit location type for levels L = 0, 1,2,3, oo. 
The level- 1 pseudothreshold is about 4.6 times greater than the asymptotic one-qubit gate 
threshold since the replacement R(l) includes many wait locations, whose initial setting is 
one-tenth of the initial one-qubit gate failure probability. The Monte-Carlo simulation gives 
a level- 1 pseudothreshold of approximately 1.3 x 10 -3 versus an estimate of 1.4 x 10~ 3 from 
the flow map. 

Similarly, for a two-qubit gate location, the level- 1 pseu dothreshold is a factor of 2 larger 
than the asymptotic two-qubit gate threshold (Figure TT2|) . Note that the level- 1 two-qubit 
gate pseudothreshold is about half the size of the level- 1 one-qubit gate pseudothreshold. 
This is because error correction is required on two logical qubits and thus there are twice 
the number of locations in a one-qubit gate replacement. The level- 1 pseudothreshold here is 
about 7.3 x 10~ 4 versus a Monte-Carlo estimate of about 6.4 x 10 -4 . 

From the threshold reliability information plots (TRIPs), it is apparent that multi-parameter 
maps and higher levels of concatenation are required to determine a threshold result. Across 
location types, the largest level- 1 pseudothreshold is approximately 40 times larger than the 
smallest asymptotic gate threshold. The smallest asymptotic gate threshold is the memory 
threshold, so it is appropriate to scale this gate threshold by 10 to eliminate artifacts from 
the setting. The pseudothreshold-threshold factor then becomes 4.6 for this example. 

4.3 TIFDs for the [[7, 1,3]] code 

We use the TIFD to characterize the flow of the maps based on the semi-analytical methods 
of 0. By using a TIFD instead of the TRIP, the flow of the failure probabilities as well 
as pseudothresholds can be visualized. In the quantum case, however, the TIFD is a 4- 
dimensional flow that is challenging to visualize. Instead, we take 2-dimensional projections 
to determi ne the f low. 

Figures H~3Tll4l show TIFDs involving location types I — 1,2,10. Figure H~3l shows the 
vector field r 1 ('f) — 7 projected onto the 71 -7™ plane. Note that the flows are partitioned by 



K. M. Svore, A. W. Cross, I. L. Chuang, and A. V. Aho 




Fig. 11. TRIP for a one-qubit gate location for L = 0,1,2,3, oo for the initial setting 7" = 
1/10(7?), where i indicates all location types except for a wait location. Circles indicate pseu- 
dothresholds and an asterisk marks the threshold for this gate and setting. The pscudothresholds 
occur at probabilities 1.4 X 10~ 3 , 6 X 10~ 4 , and 5 X 10~ 4 . The one-qubit gate threshold for this 
setting is 3 x 10~ 4 . 




Fig. 12. TRIP for a two-qubit gate location for L = 0, 1, 2, 3, 00 for the initial setting 7JJ, = 
1/10(7?), where i indicates all location types except for a wait location. Circles indicate pseu- 
dothresholds and an asterisk marks the threshold for this gate and setting. The pscudothresholds 
occur at probabilities 7.3 X 10~ 4 , 4.9 X 10~ 4 , and 4.2 X 10~ 4 . The two-qubit gate threshold for 
this setting is 3.4 X 10~ 4 . 
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Fig. 13. TIFD projected onto the 71-7^, plane, where y? = for all I except 7^ and 7°. The 
arrows represent the vector field flow T 1 (7) — 7. The thick lines illustrate the separatrices. The 
horizontal separatrix appears to have zero slope, but it intersects the 7j-axis at about 10 — 2 . The 
asterisk marks the wait gate asymptotic threshold for the setting in which y w = 7 and 7i = for 
all other location types. 



separatrices shown by the thick black lines. The map is independent of concatenation level in 
our approximation, so this flow fully characterizes the behavior of the map. The j w threshold 
found by the horizontal separatrix appears around 1.1 x 10~ 4 . The asymptotic threshold 
for the map restricted to the 7u,-axis is indicated by an asterisk on the 7,„-axis. Note the 
separatrix indicates a wait location pseudothreshold. 

Figure ITU is the TIFD projected onto the 71-72 plane. Again, the flows form a separatrix 
around 72 = 2 x 71. This is because there are two error correction routines in a two-qubit 
gate replacement, compared to only one error correction routine in the replacement for a 
one-qubit gate. The 72 threshold appears along the other separatrix around 2.3 x 10~ 3 . The 
asymptotic threshold restricted to the 72-axis, indicated by the asterisk, is a factor of 3.8 
below the pseudothreshold. 

In the classical setting, the TIFD indicated the asymptotic threshold, since the map was 
exact and only two-dimensional. However, it is evident from the TIFDs for the [[7,1,3]] 
code that 2-dimensional projections of the flow are insufficient to determine the quantum 
fault-tolerance threshold set T. Since the threshold set is a multi-dimensional surface, the 
two-dimensional projection fails to indicate flow in the other dimensions. Although it appears 
the separatrices indicate a separation between points that flow to zero and those that flow to 
one, it cannot be used to determine the threshold, but it may be used to determine an upper 
bound on the threshold for this example. 

5 Techniques for Determining the Asymptotic Threshold 

In Section 0] low-dimensional projections of the flow using a TIFD were used to estab- 
lish pseudothresholds. However, the TIFD could not be used to determine the asymptotic 
threshold jth- It may be possible, though, to bound 7^ for a particular map by restricting 
the map T to the axes. In this section, we describe a possible technique for upper bounding 
the fault-tolerance threshold 7^. 

Consider the following setting, the axis setting, where every initial failure probability is 
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Fig. 14. TIFD projected onto the 71^72 plane, where 7 ; ° = for all I except 7J and 7J. The arrows 
represent the vector field flow r 1 (7) —7. The thick lines illustrate the separatrices. The horizontal 
separatrix appears to have zero slope, but it intersects the 7^-axis at about 10 — 2 . We do not 
show this intersection because the intersection with the 72-axis at about 2.3 X 10 — 3 has a more 
significant role in determining the threshold. The asterisk marks the two-qubit gate threshold for 
the setting in which 72 = 7 and 7^=0 for all other location types. 



0, except for the axis of interest, i.e., 5(7) = (7, 0, 0, . . . , 0), and 7 is assigned to the location 
axis of interest. We conjecture that the level-1 pseudothreshold for this setting upper bounds 
1th- 

Consider the plot shown in Figure El of the approximated 7™ pseudothresholds of the 
[[7, 1,3]] code for the axis setting. Note that the pseudothresholds for the wait location for 
the axis setting are strictly decreasing toward a threshold. We find similar behavior for the 
other location types as well. These pseudothresholds and the threshold found for the wait 
location are lower than the pseudothresholds and thresholds for the other location types in 
the axis setting. This suggests that the L = 1 wait location pseudothreshold for the axis 
setting, the smallest level-1 axis pseudothreshold, is an upper bound on 7^. 

This conclusion that pseudothresholds with the axis setting provide an upper bound on 
the fa ult-t olerance threshold is supported by numerical evaluation of the threshold set T in 
Figure E| This figure shows four convex hulls: one pseudothreshold hull and three threshold 
hulls. The pseudothreshold hull is determined by the 7*, 72, and 7^, axis pseudothresholds, 
which are plotted as circles, while the threshold hulls were determined numerically from the 
flow map for a grid of parameters. The region of parameter space above the pseudothreshold 
hull is strictly above threshold for any value of 71m • The largest threshold hull corresponds 
to 7° m = and all points beneath it are below threshold. Similarly, the other two threshold 
hulls correspond to 7° m = 1.5 x 1CT 3 and 7° m = 3.5 x 10~ 3 . 

Interestingly, as our choice of language indicates, T appears to be a convex set equal 
to the convex hull of the axes thresholds and the origin. The edge length of the largest 
cube in T is approximately jth ~ 8.8 x 10 -5 . Furthermore, T appears to be contained 
in the convex hull of the axes pseudothresholds and the origin. These pseudothresholds 
are -fl w 4.4 x 10~ 2 , 72 ~ 2.3 x 10~ 3 , and 7^ w 1.1 x 10~ 4 for the corresponding axes 
settings. The jl m pseudothreshold is comparable to 7*. The smallest level-1 pseudothreshold 
is 7^ Ri 1.1 x 10~ 4 , so this is an upper bound on 7^ for this example. Though the error 
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Fig. 15. TRIP for a wait location for L = 0, 1, 2, 3, 4, oo for the axis setting. The pseudothresholds 
are 1.1 X 10" 4 , 9.6 X 10" 5 , 9.2 X 10" 5 , and 9.0 X 10" 5 . The threshold for the wait location for 
this setting is 8.8 X 10 -5 . 




Fig. 16. The set T for the [[7, 1, 3]] code together with a slice of the convex hull of the level-1 axis 
pseudothresholds. The levcl-1 axis pseudothresholds are plotted with open circles and connected 
with lines to illustrate their convex hull. The other three hulls are schematic representations of the 
numerically computed boundary of the set T for varying values of 7P . All points beneath a given 
hull are below threshold. The largest cube contained in T has edge length "f th fa 8.8 X 10 -5 . Note 
that this does not contradict the wait location threshold in Figure ITU! because the other locations 
are 10 times less reliable in that calculation. 
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correction networks are slightly different, this upper bound does not contradict a rigorous 
lower bound of 2.73 x 10~ 5 for the same code [TT)] . 

To further confirm the results found using the semi- analytical map, we use a Monte- 
Carlo simulation to determine pseudothresholds for the [[7,1,3]] code. We find the following 
pseudothresholds with the axis setting by fitting a quadratic to the numerical TRIPS: 7* ~ 
8.0 x 10 -2 , 72 ~ 1.7 x 10~ 3 , and 7^ « 1.5 x 10~ 4 . The 7* pseudothreshold differs by a factor of 
1.81 from the 7J pseudothreshold found using the flow map. The 72 and 7^ pseudothresholds 
found using the flow map differ by a factor of 0.73 and a factor of 1.36, respectively, from 
numerical calculations based on a Monte-Carlo simulation. These differences could be reduced 
by revisiting some of the approximations in the flow map derivation. 

While we do not prove our conjecture, we offer two supporting observations. The first 
is an observation that follows from the fact that the flow map has a threshold. There are 
positive integers A and t such that r](7) < Ay^x f° r au A where jmax = max 7^. These 
integers determine the well-known lower bound A^ 1 /^ 1 ^ = 7™" < jth on the threshold. 

The next observation is that if T causes all components of 7 to increase or remain un- 
changed, then 7 is above the established lower bound 7™". More precisely, if 1^(7) > j£ 
for all £, then 7^ > 7™" for at least one I. This is true because if 7^ < 7™" for all £, then 
Imax < 7™"- In particular, Tl((j max ,j max , . . .,j max )) < Aj^ x < j max . 

A 7 satisfying the second observation is not necessarily above threshold, but we conjecture 
that this second observation remains true when 7™™ is replaced by "fth for maps T describing 
failure probabilities under independent stochastic error models. We know this to be the case 
for one-dimensional maps. If the conjecture is true in general, then the level- 1 pseudothreshold 
for the axis setting upper bounds 7^ for a particular map. 

6 Conclusions and Future Work 

Pseudothresholds are the failure probabilities below which recursive simulation improves 
the reliability of a particular component. Yet, pseudothresholds have been shown to be up 
to a factor of 4 greater than the asymptotic threshold for the Steane setting and more than 
a factor of 10 different for the axes settings. This behavior is a generic phenomenon in both 
classical and quantum fault-tolerance. The tools we have presented provide a way to visualize 
pseudothresholds and thresholds, and we conjecture that for a given setting, they may provide 
an upper bound on the fault-tolerance threshold for a particular map. 

Pseudothreshold behavior also influences the accuracy of some quantum threshold esti- 
mates. If some of the reported threshold estimates are actually pseudothresholds, the ex- 
amples we have given suggest that these estimates may be inaccurate by a factor of 4 or 
more. These observations apply, in particular, to some numerical threshold estimates. How- 
ever, other factors such as the noise model and circuit construction more greatly influence 
the threshold value and by making a judicious choice of concatenation level and setting, the 
difference between estimates and an asymptotic threshold can be reduced. 

Although a fault-tolerance threshold for infinite scalability cannot be determined by low- 
level pseudothresholds, pseudothresholds will become important design parameters in engi- 
neering a quantum computer. In practice, quantum computers may operate very close to 
threshold and require only a few levels of recursive simulation. If this is the case, then pseu- 
dothresholds can help determine design trade-offs and the required relative reliability of circuit 
components. In addition, the difference between pseudothresholds can be used to determine 
an appropriate level of concatenation that is within current physical capabilities. 

It is important to note that pseudothresholds can also be used to determine the frequency 
of error correction required for certain location types. If, for example, a level-1 pseudothresh- 
old for a wait location shows the failure probability worsens upon concatenation, then it may 
be beneficial to not error correct each wait location, or to only error correct a wait location 
at higher levels of concatenation. This scheme will still demonstrate a threshold; however the 
failure probability map will have to take into account the differences between the replacement 
rules for each type of location at different levels of concatenation. By optimizing the fre- 
quency of error correction at different levels of concatenation and for different location types, 
the fault-tolerance threshold can be improved. 
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This work can be extended in several directions. First, we have put forward a conjecture 
that pseudothresholds may lead to upper bounds on "f t h- If this conjecture is true, then it 
would be interesting to determine how the concatenated flow map formalism could be modified 
to give rigorous bounds. It would also be useful to determine how much a pseudothreshold 
can differ from the fault-tolerance threshold. 

Second, the analyses presented here only account for Clifford-group gates. We did not 
analyze a Toffoli, 7r/8, or other nontrivial gate needed for computational universality. It is 
still possible to apply the methods in this paper to those gates, but the corresponding flow 
map component is more difficult to estimate well. To determine a quantum fault-tolerance 
threshold, it is necessary to evaluate a computationally universal basis. With a universal 
basis, how much does the fault-tolerance threshold and sequence of pseudothresholds change? 
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